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Geometric quantum computation is a scheme to use non-Abelian Holonomic operations ratlier 
than the conventional dynamic operations to manipulate quantum states for quantum information 
processing. Here we propose a geometric quantum computation scheme which can be realized with 
current technology on nanoscale Josephson-junction networks, known as a promising candidate for 
solid-state quantum computer. 
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I. INTRODUCTION 

The elementary units of quantum information process- 
ing are quantum two-state systems, called quantum bits 
or "qubits". Unhke a classical bit, a qubit can be in 
any superposition a|0) -I- /3|1) (with a and (3 arbitrary 
complex numbers satisfying the normalization condition 
|ap -I- — 1) of the computational basis states |0) 
and |1). A qubit needs not only to preserve quantum 
coherence for a sufficiently long time but also to allow 
an adequate degree of controllability. Among a number 
of ideas proposed so far to realize qubits, the ones based 
on solid-state devices have attracted interest due to the 
scalabihty for massive information processing, _which can 
make a quantum computer of practical value J,] . 

Another crucial element of quantum information pro- 
cessing is the ability to perform quantum operations on 
qubits in a controllable way and with sufficient accuracy. 
In most proposed schemes such quantum operations are 
unitary, and conventionally have been achieved based on 
quantum dynamics governed by the Schrodinger equa- 
tion. 

Recently, it has been proposed that controllable quan- 
tum operations can be achieved by a novel geometric 
principle as well 0, Q- When a quantum system un- 
dergoes an adiabatic cyclic evolution, it acquires a non- 
trivial geometric operation called a holonomy. Holonomy 
is determined entirely by the geometry of the cyclic path 
in the parameter space, independent of any detail of the 
dynamics. If the eigenspace of the Hamiltonian in ques- 
tion is non-degenerate, the holonomy reduces to a simple 
phase factor, a Berry phase. Otherwise, it becomes a 
non-Abelian unitary operation, i.e., a non-trivial rota- 
tion in the eigenspace. It has been shown that universal 
quantum computation is possible by means of holonomies 
only 0, Q- Further, holonomic quantum computation 
schemes have intrinsic tolerance to certain types of com- 
putational errors 

A critical requirements for holonomies is that the 
eigenspace should be preserved throughout the adiabatic 
change of parameters, which is typically fulfilled by sym- 
metry 1^. It is non-trivial to devise a physical system 
with a proper eigenspace which will serve as a compu- 
tational space. Recently a scheme for geometric quan- 
tum computation with nuclear magnetic resonance was 



proposed and demonstrated experimentally Q- A sim- 
ilar scheme was proposed on superconducting nanocir- 
cuits 0. In these schemes, however, geometrically avail- 
able was only Abelian Berry phase, and additional dy- 
namic manipulations were required for universal quan- 
tum computation. A scheme based solely on holonomies 
has been proposed for quantum optical systems @ . How- 
ever, it relies on nonlinear optics, which may make this 
scheme less practical. More recently, another holonomic 
quantum computation scheme has been proposed for 
trapped ions ^1^]. Supposedly, it is the only experi- 
mentally feasible scheme proposed so far for holonomic 
computation. Here we propose a scheme for holonomic 
quantum computation on nanoscale Josephson networks, 
known as a promising candidate for solid-state quantum 
computer 0, 0, 0. It relies on tunable Josephson 
junction and capacitive coupling, which are already vi- 
able with current technology. 



II. JOSEPHSON CHARGE QUBITS 

A "Josephson (charge) qubit" [ill IT^ can be realized 
by a small superconducting grain (a Cooper-pair box) of 
size ~ lOOnm, coupled to a large superconducting charge 
reservoir or another Cooper-pair box via a Josephson 
junction, see Fig. ^ The computational bases are en- 
coded in two consecutive charge states |0) |0) and 
|1) |1) with \n) denoting a state with n excess Cooper 
pairs on the box. States with more (or less) Cooper pairs 
are suppressed due to the strong Coulomb repulsion (the 
gate-induced charge 2eng is tuned close to 1/2 mod 1), 
characterized by the charging energy Ec — (2e)^/2C 
(with C the total capacitance of the box). Excitation 
of quasiparticles is also ignored assuming sufficiently low 
temperature. The tunneling of Cooper pairs across the 
junction, characterized by the Josephson coupling energy 
Ej (^ Ec): allows coherent superposition of the charge 
states |0) and |1). A tunable junction is attained by two 
parallel junctions making up a SQUID (superconduct- 
ing quantum interference device) with a magnetic flux $ 
threading through the loop, see upper panel of Fig.lD 14j. 
Namely, in the two-state approximation, the Hamilto- 
nian is written in terms of the Pauli matrices and 
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FIG. 1: A schematic of a "Josephson charge qubit". The cross 
denotes a tunable Josephson junction, which consists of two 
parallel tunnel junctions (upper panel). The strength of the 
effective Josephson coupling is tuned by flux threading the 
loop. The gate voltage Vg controls the induced charge on the 
box. 

fT± = ^{a''±iay) as 

H = -h{n,)a^ - 1 [J{cb)a+ + J*(0)a-] (1) 

where h{ng) = Eci^rig — 1), J(0) is the effective Joseph- 
son coupling of the tunable junction (i.e. SQUID), and 
(j) = 7r$/$o (we assume that |0| < 7r/2) with $0 = hc/2e 
the superconducting flux quantum. Given Josephson en- 
ergies Ej and -fEj (7 > 0) of the two parallel junctions 
on a SQUID loop, the magnetic flux gives rise to a phase 
shift a(4>) as well as an amplitude modulation ^(0) of the 
effective Josephson coupling J(0) = 2Ej A{(f>)e~^°''^'^\ 
A{(j)) and «(</>) are given by 8] 

A{(t)) = V(l-7)V4 + 7Cos2 (2) 

and 

1 - 7 

tana{(f)) = -^-^ — tancj), (3) 

respectively. It is worth noticing that for identical junc- 
tions (7 = 1), (i) there is no phase modulation [«((/>) = 0] 
and (ii) the effective Josephson coupling can be turned 
off completely [J('?I>) = 0] at = 7r/2. In what follows, 
some tunable junctions have 7 = 1 and others 7 7^ 1 
depending on their roles for the system. 

Below we will demonstrate that one can obtain the 
three unitary operations, Uz{^) = exp (z(y9|l) (1|) (phase 
shift), Ux{^) = exp {iLpa^) (rotation around x axis), and 
Ucziv) — exp (-|-i(^|10)(10|) (controlled phase shift), on 
an arbitrary qubit or pair of qubits, using geometric ma- 
nipulations only. It is known that these unitary opera- 
tions form a universal set of gate operations for quan- 
tum computation 0, ^|. Since the charge degrees of 



freedom is used in the present scheme, the state prepa- 
ration and the state readout, which are another impor- 
tant procedures required for quantum computation, can 
be achieved using the same methods used in dynamical 
schemes [Tll |. 

III. ELEMENTARY GATES 

Before demonstrating the geometric implementation of 
elementary gates, we suggest a prototype Hamiltonian 
which reveals the proper symmetry for geometric manip- 
ulations in question. All the Hamiltonians considered in 
this paper share the following common structure: 

H = e\6){d\--Y,(nS){o\ + ^m(i\) ■ (4) 

Here fi.; is the transition amplitude from the state |0) to 
\i) and e is the energy of the state |0) measured from the 
degenerate energy of the states \i) {i = 1, . . . , N). For 
our consideration below, one may regard the state vector 
\i) (« = 0, 1, . . . , A^) represent, e.g., an excess Cooper pair 
on the zth superconducting box, see Fig. [21(a). \i) may 
also represent electronic levels in atoms such as discussed 
in Ref.E 

As one changes the parameters the Hamiltonian 
in Eq. Q preserves the (A^— l)-dimensional degenerate 
subspace. This can be clearly seen by deflning = 
^'~^Y.f=i^^\i) with = ^|f]i|2 + . .. + |f7^|2, and 
rewriting H as 

H = e\b){b\-hi(p){Q\ + \bm) . (5) 

The Hamiltonian in Eq. ((JJ corresponds to a particle in 
a (biased) double-well potential in the tight-binding ap- 
proximation. Therefore, it has two eigenstates |A±) = 
T |6)) with energies A± = l{e±yWT^). The 
other A^ — 1 levels out of A^ -|- 1 form a degenerate sub- 
space £n-i with energy 0, which we will use later for a 
computational subspace, see Fig. |21(b). Notice that the 
degenerate eigenspace £n-i is always perpendicular both 
to |0) and as the parameters fi^ change, \Q) rotates 
in the Hilbert space, and the eigenspace £n-i is attached 
rigidly, perpendicular to \Q), see Fig. |21(c). 

We first show how to get the unitary operation Uz geo- 
metrically. We consider three Cooper-pair boxes coupled 
in series as shown in Fig. 01 The first (solid) box en- 
codes the computational bases while the second and third 
(empty) boxes serve as "ancilla" qubits. The Hamilto- 
nian is given by [see Eq. (^] 

H^--Y, {Jn<ji^o + h.c) - -ha^ (6) 

n=l 

with h.c. meaning the Hermitian conjugate. Compar- 
ing Eq. O) with the prototype Hamiltonian in Eq. I^J, 
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FIG. 2: (a) A schematic representation of the prototype 
model, Eq. which has a proper degenerate eigenspace for 
geometric quantum computation, (b) The degenerate level 
structure of the model, (c) The structure of the correspond- 
ing Hilbert space, in which the degenerate subspace is always 
perpendicular to \Q). 



noticed are following correspondences: |6) —> 1 100) 012 
(|100)i23 is short for |l)o|0)i|0)2), |i) ~^ |010)oi2, |2) 
1001)012, |3) -> 1000)012, ^^l ^ Ji, ^2 ^ J2, and ris = 0. 
From this [or direct diagonalization of the Hamiltonian 
Eq. lO, of course], one can see that the two states 

I Ai) = J* 1010)012 -Ji* 1001)012 (7) 



(not normalized) and 



IA2) = |000)oi 



(8) 



form a degenerate subspace with energy —h/2, which is 
preserved during the change of Ji and J2 (equivalently (pi 
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FIG. 3: Configuration for the phase-shift operation Uz- 
Throughout the work, computational bases are encoded in 
the solid boxes and empty boxes serve as ancilla qubits. 



and 4>2)- Since the computational basis is only encoded 
in the "true" qubit (box 1), the total wave function |\1/) 
of the logical block should be initially prepared in a sepa- 
rable state with respect to the true qubit and the ancilla 
qubits, 1^*) = 17/1)1 ® \4'')q2- In other words, one should 
be able to turn off at will the tunable junction 1, which 
should therefore be made of identical parallel junctions 
(71 — 1), see Eqs. (0) and JSJ. After a cyclic evolution 
of the parameters (pi and 02 along a closed loop starting 
and ending at the point with (pi = 7r/2 (i.e., Ji = 0), the 
state |Ai) acquires the Berry phase [I^. For example, 
along the loop depicted in Fig.|31(b), the Berry phase is 
given by 
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2) 



1^ 



9) 



(The Berry phase vanishes if 72 = 1 as expected 1(^.) 
The state IA2) remains unchanged. Therefore, the cyclic 
evolution amounts to Uz{'fB)- 

Here it should be emphasized that although used is 
the Abelian Berry phase, the degenerate structure is cru- 
cial. The dynamic phases acquired by |Ai) and IA2) are 
the same and result in a trivial global phase. In recently 
proposed schemes 0, |^ , which have no degenerate struc- 
ture, at least one dynamic manipulation was unavoidable 
to remove the dynamically accumulated phases. Another 
point to be stressed is that the phase shift a2{4>2) in the 
effective Josephson coupling J2{4>2) is indispensable for 
the Berry phase. 

The two-qubit gate operation Ucz can be realized 
geometrically using capacitive coupling. As shown in 
Fig- El (a), the ancilla qubits on different three-box sys- 
tems are coupled in parallel via capacitors with capaci- 
tance Cc- It is known [H El El that for Cc suf- 
ficiently larger than the self-capacitance C of each box, 
the states |01)oo' and |10)oo' are strongly favored over the 
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FIG. 4: (a) Configuration for the two-qubit gate Ucz- (b) Ef- 
fective "joint-tunnel" coupling (crossed wiggly lines) resulting 
from the capacitive coupling in (a). 



states |00)oo' and |ll)oo') and the same for boxes 2 and 
2' (recall that Ug ^ 1/2 for each box). This effectively 
leads to a "joint tunneling", see Fig. ^ (b): Tunneling 
of a charge from box 1 to should be accompanied by 
tunneling of another charge from 0' to 1'. The joint- 
tunneling amplitudes are given by Ja ~ AJiJ*,/Ec and 
Jb ~ AJ2J21/EC 0. Then the total Hamiltonian has 
the form 




FIG. 5: (a) Configuration for a rotation around x axis, Ux- 
The computational bases is encoded over the two (solid) 
boxes, not in a single box. (b) An adiabatic path in the 
parameter space to achieve Ux {4>3 is kept constant all along 
the path). 



- I [Jb (atcT^) (a2;4) + h.c] - [a'o - <'] ■ (10) 

In analogy with Eqs. Q and ®, the above Hamiltonian 
has an eigenspace containing the four degenerate states 

|Aoo) = |00)n' ^ |0101)oo'22' , (11) 
|Aoi) = |01)n' ® |0101)oo'22' , (12) 
|Aii) = |ll)ir®|0101)oo'22', (13) 

and 

|Aio) = Jfc |10)ir |0101)oo'22' - J:|01)ii' ® |0110)oo'22' 

(14) 

(not normalized) with energy ~h. As in the previous case 
[see discussions below Eq. (jHl], it is assumed that the 
tunable junctions Ji and Ji' are composed of identical 
parallel junctions (71 =71' = 1), while 72 ^ 1 and 72' ^ 
1. One can fix (jjii — 02' = and change the parameters 
(pi and 4>2 along a closed loop starting from the point with 
01 = 7r/2 (Ji = Ja = 0). Upon this cyclic evolution, the 
state |Aio) acquires the Berry phase in Eq. ^ while the 
other three states remain unchanged, leading to the two- 
qubit gate operation Uczifs)- We mention that in this 



implementation of Ucz, the capacitive coupling is merely 
an example and can be replaced by any other coupling 
that effectively results in a sufficiently strong Ising-type 
interaction of the form Jija^a^. 

In the demonstrations of implementing Uz and Ucz 
above, we have encoded the bases |0) and |i) in a single 
Cooper-pair box for simplicity. To realize Ux , we need to 
encode the basis states over two Cooper-pair boxes, e.g., 
box 1 and 2 in Fig.El(a): |0) |01)i2 and |I) |10)i2. 
It is straightforward to generalize the above implemen- 
tations of Uz and Ucz in this two-box encoding scheme. 
Now we turn to the remaining single-qubit gate operation 
Ux- The Hamiltonian is given by 

H=-l {Jn<ji<Jo +h-c.)-ha^. (15) 

n=l,2,3 

The degenerate subspace is defined by the two eigenstates 

|Ai) = (J2* |10)i2 - Ji 101)12) ® |00)o3 (16) 
(not normalized) and 

1^2) = I , J p (Jl|10)l2 + J2|01)l2) ® |00)03 

\Ji\ +\J2r 

-|00)l2<»|01>03 (17) 

(not normalized) both with energy — /i/2. In this case, it 
is required that 71 = 72 = 1 but 73 7^ 1 [see discussions 
below Eq. ©]. As an example, we take a closed loop 
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FIG. 6: A computational network for universal holonomic 
quantum computation. 



shown in Fig. El (b) (one may choose any path starting 
and ending at 4>i = 4*2 = 7r/2, i.e., Ji = J2 = 0) with (/)3 
fixed. The adiabatic theory of holonomies Q ensures that 
from this adiabatic cycle, a state \tp) initially belonging to 
the eigenspace undergoes a change to \tp') = Ultp). The 



unitary operator U is given by U 
with (fi' — a3((/)3)/2 — 7r/4 and 



Ul{ip')Ux{v)Uziip') 



dx cos ( 



0. (a;2 + cos2 + (a;2 + cos^ (?!),)Mi(03) 

(18) 

Removing the first and last factors oi Uz (if necessary) 
with additional phase-shift operations, one can achieve 
Ux- 

Finally, a quantum network can be constructed as in 
Fig. |S1 to perform all the unitary operations discussed 
above (and hence universal quantum computation) geo- 
metrically. It is noted that the coupling capacitance Cc 
is not tunable, but it suffices to have a control over each 
tunable junction (i.e. SQUID) and gate voltage. 



IV. DISCUSSION 

So far we have discussed fundamental requirements for 
geometric manipulations in idealized systems. In this 
section, we discuss several situations one may come up 
with when attempting an experimental realization of the 
present scheme. 

First of all, t he p resent scheme is based on the adi- 
abatic theorem zl]. Ideally the change in the control 
parameters should be infinitely slow. At a finite rate 
of change, there can be a transition out of the compu- 
tational subspace. However, typically such a Landau- 
Zener-type transition occurs with an exponentially small 
probability P ~ exp(— ttA/t;), where ?iA is the smallest 
energy gap between the computational subspace and the 
nearby energy levels and 7 7 is the adiabaticity parameter 
(i.e., §iH{t) ~ Ty/iA) [H lH H 11 11 • In Josephson 
networks the level distance is of order of the Josephson 



coupling energy, A Ej/h ~ (80ps)~^. For an opera- 
tion time Top = l/?7 ^ 3/ A, P < 10^'*. (For compari- 
son, in a recent experiment concerning dynamic quantum 
computation on Josephson qubits , the switching time 
was no shorter than 0.2511/ Ej ~ 20 ps.) 

In the previous section, on each superconducting box 
we neglected higher charge states other than the two low- 
est. In dynamic quantum computation schemes, the exis- 
tence of those higher levels may cause the quantum leak- 
age errors, i.e., it leads to a nonzero probability of leak- 
age out of the computational space, and more severely to 
renormalization of the energy levels in the computational 
space (which therefore reduces the gate fidelity) In 
Josephson qubits, however, the coupling to the higher 
charge states are only through the Josephson tunneling 
of Cooper pairs, which can be easily included in our con- 
siderations [see Eqs. Q, 10), (jTHI), and ^^]. Those 
higher charge states form well-separated energy levels, 
and do not alter the degenerate structure of the sub- 
space in question, at least up to the order of {Ej/EcY' 
[in the dynamic scheme the quantum leakage error oc- 
curs already in the order of {Ej/Ec}'^, see Ref. 27J. The 
leakage to the higher charge states out of the compu- 
tational subspace can therefore be considered within the 
framework of Landau-Zener tunneling, which has already 
been discussed above. The renormalization of (degener- 
ate) energy of the computational space is irrelevant in 
our geometric scheme since it does not rely on the dy- 
namical time-evolution operator but only on the purely 
geometric means. 

In reality there are fluctuations of the (reduced) flux cj) 
(tuning the junctions) and the gate-induced charge 2eng 
(resulting from the fluctuations of random charges in the 
substrate or gate voltage itself). One consequence of 
these fluctuations is the Landau-Zener-type transitions 
out of the computational subspace. A recent experiment 
on Josephson charge qubits [sj suggests that fluctua- 
tions of Ug as well as (f) are dominated by low frequency 
fluctuations. Therefore, the Landau-Zener-type transi- 
tions might be small. The fluctuations of rig can cause 
another type of errors: While the eigenspace is by con- 
struction robust against the low-frequency fluctuations 
of 0, the random charge fluctuations lift the degeneracy 
of the computational subspace. The wave function of the 
system then acquires dynamically accumulated phase fac- 
tors ex-p{—iSETop/h), where SE is the small level spac- 
ing caused by the fluctuations of Ug. Such dynamical 
phases can be ignored for sufficiently small fluctuations 
and sufficiently short - yet long enough for adiabaticity 
- operation time (1/A ^ Top <C 1/SE). 

Another common source of decoherence in Josephson 
charge qubits is the quasi-particle tunneling '2^ . In par- 
ticular, since the computational eigenspace is not the low- 
est energy states [Eq. (0J and Fig. |21 (b)], it gives rise 
to the relaxation out of the eigenspace to lower energy 
states (this effect cannot be described by the Landau- 
Zener-type transitions). At sufficiently low temperatures 
compared with the superconducting gap A5, the quasi- 
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particle tunneling rate Fqp is exponentially small [2^ . 
Fqp ~ exp[— (2A5 + Ec)/kBT]. For example, in the 
experiment on Cooper-pair box Fqp ^ (6ns)~^ ~ 
10~^Ej/h at temperature 30 mK even through the probe 
junction which was biased by a voltage eV ^ 2A5 + Ec 
(without the voltage bias Fqp should be even smaller). 
In such a situation, one can therefore conclude that the 
effect of quasi-particles is negligible. 

For a brief comparison of the present scheme with the 
conventional Josephson charge qubit we estimate 

the fidelity for a single phase-shift operation Uz- The 
fidelity in this case is given by 

fidelity ~ y^(l - P)[l - sin'*(J0/2)] , (19) 

where P is the probability for Lanau-Zener-type tran- 
sitions or quasiparticle tunnelings to occur and (50 is 
the error in phase shift due to the background charge 
fluctuations, i.e., Sip = SETop/h (see above). Taking 
A = (80ps)-i, Top = 3/A, As = 5A, fc^T = A/10, 
and 6E/T1 = A/10 (see Ref. 12), one estimates fidelity ~ 
0.998. In the dynamic scheme the fidelity is also given by 
the same form as Eq. (|19|) . The differences are that P is 
mainly responsible for the quasiparticle tunnelings and 
that the phase error S4> comes from the finite ramping 
time of the gate pulse. Taking the parameters from the 
recent experiment |0| , we see that the fidelity takes the 
same value (up to three digits below decimal point). 

Lastly, in ideal case some tunable junctions [e.g., Ji in 
Eq. ©, see discussions below Eq. ©] need to be turned 
off completely. In reality, the Josephson energies of the 
two parallel junctions on a SQUID loop (Fig.^l may not 
be identical [i.e., 7 7^ 1 in Eq. ^ and Fig. [J. Then, 
a tunable junction (i.e. SQUID) cannot be turned off 
completely 0. This makes it nontrivial to prepare an 
initial state which should be a product state of the "true" 
qubit and the ancilla qubits in a logical block [see, e.g., 
|Ai) below Eq. ©]. In practice, such a difficulty can 
be overcome by means of fast relaxation processes with 
the gate voltages of the ancilla qubits adjusted far off 



the resonance in the initial state-preparation stage. This 
process also allows for preparation of the "true" qubit in 
a definite initial state llll. 



V. CONCLUSION 

We have proposed a scheme based on geometric means 
to implement quantum computation on solid-state de- 
vices. The main advantage of a geometric computation 
scheme is its intrinsically fault-tolerant feature 0, 01 ■ 
However, it is usually non-trivial to find a physical system 
whose Hamiltonian has a particular degenerate structure 
for geometric computation. The scheme discussed in this 
work provides a generic way to construct such a system 
from arbitrary quantum two-state systems as long as cou- 
plings satisfy certain requirements specified above. Such 
requirements are rather easy to fulfill on solid-state de- 
vices. A drawback of this scheme is that it requires more 
resources (four Cooper-pair boxes per each qubit). Con- 
sidering quantum error correcting codes, however, it may 
not be a major disadvantage. Moreover, since the cur- 
rent scheme is based on adiabatic evolution, it does not 
require sharp pulses of flux and gate voltage. With cur- 
rent technology, it is still challenging to obtain a suffi- 
ciently sharp pulses of flux and gate voltages (in Ref. 0, 
the raising and falling times were of order of h/Ej). Fi- 
nite raising and falling times of pulses can result in a 
signiflcant error in dynamic computation schemes . 
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